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Abstract 

A field-theoretic description of the critical behaviour of the weakly disordered 
systems is given. Directly, for three- and two-dimensional systems a renormal- 
ization analysis of the effective Hamiltonian of model with replica symmetry 
breaking (RSB) potentials is carried out in the two-loop approximation. For 
case with 1-step RSB the fixed points (FP's) corresponding to stability of the 
various types of critical behaviour are identified with the use of the Pade-Borel 
summation technique. Analysis of FP's has shown a stability of the critical 
behaviour of the weakly disordered systems with respect to RSB effects and 
realization of former scenario of disorder influence on critical behaviour. 



64.60.Ak, 64.60.-i, 64.60.Fr, 64.60.Ht 



1 



Typeset using REVTgX 



The effects produced by weak quenched disorder on the critical phenomena have been 
studied since many years ago According to the Harris criterion Q, the disorder effects 

the critical behaviour only if a, the specific heat exponent of the pure system, is positive. 
In this case a new universal critical behaviour, with new critical exponents, was established. 
In contrast, when a < 0, the disorder appears to be irrelevant for the critical behaviour. 

In dealing with the weak quenched disorder the traditional approach is the replica method 
and in terms of replicas all the results obtained for the disorder systems correspond 
to the so-called replica-symmetric (RS) solutions. Physically it means that only unique 
ground state is assumed to be relevant for the observable thermodynamics. However, in a 
number of papers |]-|8|] the ideas about replica symmetry breaking (RSB) in the systems 
with quenched disorder were presented. For the first time in the paper || physical arguments 
showing that in the presence of the quenched disorder there exist numerous local minimal 
energy configurations separated by finite barriers and demonstrating how the summation 
over these local minima configurations can provide additional RSB interaction potentials 
for fluctuating fields have been offered. The renormalization group (RG) description of the 
classical 4 model with RSB potentials was presented in the one-loop approximation using 
e - expansion It was shown that the RSB degrees of freedom produce dramatic effect 

on the asymptotic behaviour of the RG flows, such that for a general type of the RSB there 
exist no stable fixed points, and RG equations arrive into the strong coupling regime. In 
contrast, in Ref. |J using of the fermion representation it was demonstraited that the critical 
behaviour of the 2D weakly disordered Ising system is stable with respect to the RSB modes. 

However, our numerous investigations of pure and disordered systems performed in the 
two-loop and higher orders of the approximation for the 3D system directly together with 
methods of series summation show that the predictions made in the lowest order of the 
approximation, especially on the basis of the e - expansion, can differ strongly from the real 
critical behaviour flO|| . Therefore, the results of RSB effects investigation in Refs. |6|-§ must 
be reconsidered with the use of a more accurate field-theory approach in the higher orders 
of the approximation. 

In this paper, we have realized the field-theoretical RG description in the two-loop ap- 
proximation of the 3D and 2D models of the weakly disordered systems with RSB interaction 
potentials of forth order on fluctuating fields. We have carried out the solution of the RG 
equations with the use of series summation method and the analysis of stability of various 
types of critical behaviour with respect to RSB effects. 

We consider an 0(p) - symmetric Ginzburg-Landau- Wilson model of a spin system with 
weak quenched disorder near critical point given by the Hamiltonian 

H = J ( ltm(x)} 2 + hr - *T(aO]E#(*) + \g E , (!) 

[ 8=1 - 8=1 i,j = l ) 

where <f>%(x) is the p-component order parameter and Sr(x) is the Gaussian-distributed ran- 
dom transition temperature with the second moment of distribution (((5t(x)) 2 )) ~ u defined 
by the positive constant u which is proportional to the concentration of defects. The use of 
the standard replica trick gives the possibility to average easily over the disorder and reduce 
the task of statistical description of the weakly disordered system with the Hamiltonian ([!]) 
to the homogeneous system with the effective Hamiltonian 
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ff » = / ^ { n E E [[v#(^)] 2 + r[c(^)] 2 ] + 1 E E ^[0- (^)] 2 [05(^)] 



=1 a=l 



(2) 



i,J = l a, 6=1 



which is a functional of n replications of the original order parameter with an additional 
vertex u in the replica symmetric matrix g ab = g5 ab — u. The properties of the original 
disordered system are obtained in the replica number limit n — > 0. The following standard 
RG procedure based on the statistical calculation of contribution to the partition function of 
longwavelength order parameter fluctuations around the global minimum state with <p(x) = 
gives the possibility to find the various types of critical behaviour and conditions of their 
stability and carry out the calculation of critical exponents. 

However, as it was shown in Refs. @-f| , the fluctuations of random transition temperature 
St(x) for [t — 5t(x)\ < can lead to realization in system of numerous regions with <p{x) ^ 
displaying through the numerous local minimal energy configurations separated from the 
ground state by finite barriers. In this case the direct application of the traditional replica- 
symmetric RG scheme may be questioned. For statistical description of such systems near 
phase transition point the Parisi RSB scheme was used in Refs. [0-fl by analogy with 
spin glasses 0|. It was argued that spontaneous RSB can occur due to the interaction 



of the fluctuating fields with the local non-perturbative degrees of freedom coming from 
the multiple local minima solutions of the mean-field equations. It was shown that the 
summation over these solutions in the replica partition function can provide the additional 
non-trivial RSB potential Y^ab 9ab<f )2< f )2 in which the matrix g ab has the Parisi RSB structure 
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According to the technique of the Parisi RSB algebra, in the limit n — > the matrix 
g a b is parametrized in terms of its diagonal elements g and the off-diagonal function g(x) 
defined in the interval < x < 1: g ab — > (g,g(x)). The operations with the matrices g ab are 
given by the following rules: 



9ab 



(g k ;g k (x)),(g\ b 



J2dac9cb -> (c; c(x)), (g 3 ) ab = J2 9ac9cd9db -> (d; d(x)), (3) 



cA=l 



where 



c = g 2 — / dxg 2 (x), c(x) = 2 
Jo 



9-1 dyg(y) 
Jo 



g(x) - / dy[g(x) - g{y)]\ (4) 
Jo 



d = eg — / dxc(x)g(x), d(x) 
Jo 



9- / dyg(y) 



c(x) + 



<■- I dyc(y) 
o 



9{ x ) 



- / dy[g{x) - g{y)][c(x) - c(y)]. (5) 
Jo 



The RS situation corresponds to the case g(x) = const independent of x. 

We carried out the field-theoretical RG description of the 3D and 2D models with the 
effective replicated Hamiltonian (0) in which the matrix g ab has the RSB structure in the 
two-loop approximation. In the field-theoretic approach the asymptotic critical behavior 
of systems in the fluctuation region are determined by the Callan-Symanzik RG equation 
for the vertex parts of the irreducible Green's functions. To calculate the j3 functions as 
functions of the renormalized elements of the matrix g ab appearing in the RG equation, 
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we used the method based on the Feynmann diagram technique and the renormalization 
procedure [12|. We obtained the next expressions for two-pont vertex function and 
four-point vertex functions T^: 



dk 2 



= l+4fg 2 aa + 2pfY / 9acg ca , (6) 

fc 2 =0 c=l 



F (4) 
1 ab 



kl =o = 9ab~pJ2 9ac9cb - 4:g aa 9ab - ^gl b + (8 + 16%„ 6 + (24 + 8h)g 2 aa g ab 

c=l 



+48hg aa g 2 ab + Ag aa g bb g ab + 8ph ^ g ac g 2 ch + 8phg ab ^ 

9ac9cb 

c=l c=l 
n n n n 

+Aphg ab 9 2 ac + 2 pJ2 9ac9ccgcb + 4pg aa ^ g ac 9cb + P 2 9ac9cd9db, (7) 



c=l c=l c=l c,d=l 



where 



1 d f d d hd d k 2 

f(d) = 



J 2 dk 2 J [k 2 + l)(k 2 + l)({ki + k 2 + k) 2 + 1] 

du Ml 



k 2 =0 



J_ f d d hd d k 2 

l j W (A; 2 + l)(A; 2 2 + l)((A; 1 + fc 2 ) 2 + l)' W 
J = | A/(A; 2 + l) 2 , f(d = 3) = A ft( d = 3) = | 
/(<£ = 2) = 0.11464, = 2) = 0.78129, 

and was made the redefinition g ab — > g ab / J. 

However, the renormalization procedure for vertex functions is made difficult because 
of complicated expressions (0)-(||) for the operations with the matrices g ab . The step-like 
structure of the function g(x) revealed in Refs. ||-|| gives the possibility to realize the 
renormalization procedure. In this paper we considered only the matrices g ab which have 
the structure known as the 1 step RSB with function g(x) of the next view: 

where < xq < 1 is the coordinate of the step and it remains arbitrary parameter. The 
value of Xq does not be changed during the renormalization procedure and remains the same 
as in the starting function go(x). In consequence the RG transformations of the effective 
replicated Hamiltonian with RSB potentials are determined by the three parameters g, go 
and gi. 

The critical properties of model can be extracted from the coefficients Pi(g,go,gi) (i = 
1, 2, 3), 7,/, (p, go, gi), and 7^2 (g, g Q , g\) of the Callan-Symanzik RG equation. We obtained the 
(3 functions in the two-loop approximation in the form of the expansion series in renormalized 
parameters g, go and gi 

Pi = -9 + (8+ p) g 2 - px go 2 -p(l- x ) 9l 2 + ((8 / - 40 h + 20)p 
+ 16 / - 176 h + 88)g 3 + (24 h - 8 f - 12) x pgg 2 + (24 h - 8 / 



4 



- 12) (1 - x ) pggi 2 + (16h- 8) x pg 3 + (16h- 8) (1 - x ) Pgi 3 , 
fa = -go + (4 + 2p) gg - (2px - 4) g 2 - 2 (1 - x )pg g 1 

+ ((8 / - 48 h + 28)p + 16 / - 48 h + 2A)g 2 g Q + (((32 h - I6)x (10) 
+ 8 - 32 h)p + 48 - 96 /i)M) 2 + (32 h - 16) (1 - x ) p~gg ogi 
+ ((48 - 8 / - 20) x p - 32 h + lQ)g 3 + (32h- 8) (1 - x ) pg 2 gi 
+ (lQh-l2-8f)(l-x )pg gi 2 i 
P 3 = -g 1 -px g 2 + (p{x -2)+4) < ?1 2 + (4 + 2p)^ 1 + ((8/-48/i 
+ 28)p + 16 / - 48 h + 2A) gi g 2 + {lQh- 8) x pgg 2 + (((8 - 16 h) x 

- 8)p + 48 - 96 /i) ? i ^i 2 + (16 h - 8) x pg 3 + (8h-8f-4)x pgm 2 
+ ((8 / - 24 h + 12) x p + (48 h - 8 f - 20) p + 16 - 32 /i) ^ 3 . 



By analogy with papers ||-||[] we changed (7 a ^5 — > —g a ^b in the expressions ( |10|) for the f3 
functions, so that the off-diagonal elements g a ^b would be positively defined. 

It is well known that perturbation series are asymptotic series, and that the vertices 
describing the interaction of the order parameter fluctuations in the fluctuating region r — ► 



are large enough so that expressions ( |i~0"D cannot be used directly. For this reason, to extract 
the required physical information from the obtained expressions, we employed the Pade- 
Borel approximation of the summation of asymptotic series extended to the multiparameter 
case. The direct and inverse Borel transformations for the multiparameter case have the 
form 

oo 

f(g, go, gi) = c mfgog\ = J e^F(gt, g t, git)dt, (n) 

^^^) = E (i+ 7+ A;) , g < ^- (12) 

A series in the auxiliary variable 9 is introduced for analytical continuation of the Borel 
transform of the function: 

oo k k—i 

F(g, go, gi, 0) = £0*£E ^r^g'gy; ' '■ (13) 

fc=0 i=0 j=0 K - 

to which the [L/M] Pade approximation is applied at the point 9 = 1. To perform the 
analytical continuation, the Pade approximant of [L/l] type may be used which is known to 



provide rather good results for various Landau- Wilson models (see, e.g. Refs. | 13| , |I1 |). The 



property of preserving the symmetry of a system during application of the Pade approxima- 



tion by the 9 method, as in ||13fl , has become important for multivertices models. We used 
the [2/1] approximant to calculate the f3 functions in the two- loop approximation. 

The nature of the critical behavior is determined by the existence of a stable FP satisfying 
the system of equations 

M*,9o,9l)=Q (« = 1, 2, 3) (14) 

for resummed f3 functions. We have found three types of non-trivial FP's in the physical 
region of parameter space g*, g^, g* > for different values of p = 1, 2, 3, which are presented 
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in Tables |I]-|V| (the exception was made in the case with p = 3, when presented in the 
Table |T| the coordinates of type II and type III FP's are characterized by unphysical negative 
values of and g\). Type I with g* ^ 0, g$ = g\ = corresponds to the RS FP of a pure 
system, type II with g* ^ 0, g$ = g* ^ is a disorder-induced RS FP, and type III with 
9* 7^ 0; 9o = 0) 9i corresponds to the 1-step RSB FP. The values of parameters g*,g\ 
for the 1-step RSB FP depend on the coordinate of the step xo, and we presented in Tables 
P^jTV| the received values of these parameters in the range < xq < 1 with changes through 
the step Aa;o = 0.1. 

The type of critical behavior of this disordered system for each value of p is determined 
by the stability of the corresponding FP. The requirement that the FP be stable reduces to 
the condition that the eigenvalues Aj of the matrix 

B dm,9i9i) (15) 

dg 3 

lie in the right-hand side complex half-plane. Analysis of the received values A, for FP's in 
Tables [I-1V shows that for 3D and 2D Ising models (p = 1) and 3D XY model (p = 2) the 
disorder-induced RS FP's are stable. However, we believe that in the higher field-theory 
orders of approximation the RS FP of a pure system will be stable for the 3D XY model. 
Two facts indicate this, such as the weak stability of the disorder-induced RS FP and that 
in the two- loop approximation the marginal value of p c = 2.0114 for the borderline between 
regions of stability for the disorder-induced RS FP and the RS FP of a pure system. In 
the higher orders of approximation the marginal value of p c < 2 such as the specific heat 
exponent a > for the pure XY model. For the 3D Heisenberg model (p = 3) the RS FP 
of a pure system is stable and both another types of FP's are characterized by unphysical 
values of coordinates <?q and g\. Our conclusions coincide with results of Ref. || for the 2D 
random Ising model, but contradict with results of Refs. for 3D disordered systems. 

However, we must note that the obtained RS FP values for vertices g, go and the eigen- 
values Ai and A2 of the stability matrix correspond to results of Ref. []15| |, in which a field- 
theoretic treatment of disordered 3D and 2D spin systems was presented in the two-loop 
approximation. The vertices v 1 and v 2 introduced in Ref. [pjl are connected with the vertices 
g and go by the relations V\ — (p + 8)(g + go) and t>2 = 8go- We have calculated the static 
critical exponents from the resummed by the generalized Pade-Borel method 7 functions in 
the corresponding stable RS FP's (Table |Vp. For comparison we presented also in Table [V] 
a values of the critical exponents from [fT6|JT~7|1 received for pure and disordered 3D systems 
without RSB in the six-loop approximation. Comparison of the exponent values shows that 
their differences are not more than 0.01. It gives the possibility to consider our results of 
the RSB effects investigation as reliable. The model with RSB potentials is the another 
one example of the multivertices models [13] for which the predictions made on the basis 
of the e - expansion can differ strongly from results of the use a more accurate field-theory 
approach for the 3D system directly together with methods of series summation. This sit- 
uation is explained by a competition of numerous types of FP's in the parameter space of 
the multivertices models. Therefore, the spread of results of the e - expansion from e < 1 
to e — 1 is impossible, as a rule, without intersection of the stability ranges for the various 
types of FP's. 

Thus, the RG investigations carried out in the two-loop approximation show the stabil- 
ity of the critical behaviour of weakly disordered systems with respect to the RSB effects. 
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In dilute Ising-like systems the disorder- induced critical behaviour is realized with RS FP. 
The weak disorder is irrelevant for the critical behaviour of systems with multicomponent 
order parameter although the proof of it for 3D systems with two-component order param- 
eter demands a calculations in the higher orders of approximation. The possible influence 
the RSB degrees of freedom on the critical behaviour of highly disordered systems can be 
nonperturbatively revealed by the use of the Monte Carlo simulation method [18|] for defini- 
tion of the probability distributions for order parameter and random transition temperature 
fluctuations. 
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TABLES 

TABLE I. Coordinates of the FP's and eigenvalues of the stability matrix for the 3D Ising 



model. 



lype 




9 


9n 


9^ 


M 


^2 


A 3 


1 




0.1774103 








0.65355 


- 0.16924 


-0.16924 


2 




0.1843726 


0.0812240 


0.0812240 


0.525 ± 


0.089i 


0.211 


3 


0.0 


0.1843726 





0.0812240 


0.525319 ± 


0.089273i 


-0.039167 




0.1 


0.1839722 





0.0829404 


0.535185 ± 


0.098291i 


-0.049185 




0.2 


0.1835134 





0.0846432 


0.547065 ± 


0.106665i 


-0.059851 




0.3 


0.1829917 





0.0863186 


0.560666 ± 


0. 113305i 


-0.071187 




0.4 


0.1824035 





0.0879503 


0.576473 ± 


0. 118038i 


-0.083210 




0.5 


0.1817458 





0.0895200 


0.595060 ± 


0.120271* 


-0.095927 




0.6 


0.1810165 





0.0910067 


0.617241 ± 


0.118872i 


-0.109334 




0.7 


0.1802154 





0.0923872 


0.643936 ± 


0.111389i 


-0.123415 




0.8 


0.1793442 





0.0936384 


0.675972 ± 


0.092079i 


-0.138133 




0.9 


0.1784070 





0.0947426 


0.713456 ± 


0.035266i 


-0.153431 




1.0 


0.1774103 





0.0956920 


0.857325 


0.65355 


-0.169237 


TABLE II. 


Coordinates of the FP's and 


eigenvalues of the stability matrix for the 3D XY 


model. 
















Type 


xo 


9* 


9a 


9* 


Ai 


A 2 


A3 


1 




0.1558303 








0.667315 


-0.001672 


-0.001672 


2 




0.1558310 


0.0005837 


0.0005837 


0.667312 


0.001682 


0.000004 


3 


0.0 


0.1558310 





0.0005837 


0.667313 


0.001683 


-0.000001 




0.1 


0.1558310 





0.0006143 


0.667313 


0.001684 


-0.000088 




0.2 


0.1558310 





0.0006483 


0.667313 


0.001685 


-0.000186 




0.3 


0.1558310 





0.0006863 


0.667313 


0.001686 


-0.000296 




0.4 


0.1558310 





0.0007291 


0.667313 


0.001687 


-0.000419 




0.5 


0.1558310 





0.0007775 


0.667313 


0.001687 


-0.000559 




0.6 


0.1558309 





0.0008327 


0.667313 


0.001688 


-0.000717 




0.7 


0.1558308 





0.0008964 


0.667314 


0.001690 


-0.000901 




0.8 


0.1558307 





0.0009707 


0.667314 


0.001692 


-0.001116 




0.9 


0.1558306 





0.0010583 


0.667315 


0.001694 


-0.001369 




1.0 


0.1558303 





0.0011633 


0.667316 


0.001696 


-0.001672 


TABLE III. 


Coordinates of the FP's and eigenvalues of the stability matrix for the 3D Heisen- 


berg model. 














Type 


xo 


9* 




9*1 


Ai 


A 2 


A3 


1 




0.1382700 








0.681378 


0.131537 


0.131537 


2 




0.1419323 


-0.0358629 


-0.0358629 


0.672676 


-0.089135 


-0.005783 


3 


0.0 


0.1419323 





-0.0358629 


0.672676 


-0.089135 


-0.005783 




0.1 


0.1419931 





-0.0381865 


0.672729 


-0.086515 


0.001104 




0.2 


0.1420386 





-0.0408334 


0.672845 


-0.083560 


0.008802 




0.3 


0.1420600 





-0.0438761 


0.673046 


-0.080206 


0.017469 




0.4 


0.1420441 





-0.0474104 


0.673361 


-0.076366 


0.027299 




0.5 


0.1419699 





-0.0515650 


0.673831 


-0.071931 


0.038540 




0.6 


0.1418040 





-0.0565177 


0.674509 


-0.066754 


0.051518 




0.7 


0.1414913 





-0.0625193 


0.675476 


-0.060636 


0.066656 




0.8 


0.1409374 





-0.0699349 


0.676837 


-0.053299 


0.084517 




0.9 


0.1399720 





-0.0793131 


0.678742 


-0.044336 


0.105835 




1.0 


0.1382700 





-0.0915089 


0.681378 


-0.033124 


0.131537 
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TABLE IV. Coordinates of the FP's and eigenvalues of the stability matrix for the 2D Ising 



model. 



Type 


20 


9* 


9a 


9l 


Ai 


A 2 


A3 


1 




0.2693889 








0.461180 


-0.095180 


-0.095180 


2 




0.2735768 


-0.0645918 


-0.0645918 


0.446348 


0.162092 


0.056959 


3 


0.0 


0.2735768 





-0.0645918 


0.446348 


0.162092 


-0.025112 




0.1 


0.2733970 





-0.0665073 


0.446636 


0.165301 


-0.029995 




0.2 


0.2731825 





-0.0685156 


0.447048 


0.168803 


-0.035240 




0.3 


0.2729279 





-0.0706204 


0.447615 


0.172637 


-0.040881 




0.4 


0.2726271 





-0.0728248 


0.448377 


0.176842 


-0.046958 




0.5 


0.2722731 





-0.0751310 


0.449385 


0.181468 


-0.053513 




0.(5 


0.2718579 





-0.0775405 


0.450704 


0.186566 


-0.060592 




0.7 


0.2713724 





-0.0800528 


0.452416 


0.192194 


-0.068245 




0.8 


0.2708067 





-0.0826660 


0.454631 


0.198417 


-0.076523 




0.9 


0.2701495 





-0.0853752 


0.457490 


0.205298 


-0.085483 




1.0 


0.2693890 





-0.0881727 


0.461180 


0.212903 


-0.095180 



TABLE V. Critical exponents of the 3D models for RS FP's 



Model 


FP 




V 


V 


7 


P 


a 


Ising 


RSI 


this work 


0.028 


0.631 


1.244 


0.324 


0.107 






Ref. (l| 


0.031(4) 


0.630(2) 


1.241(2) 


0.325(2) 


0.110(5) 




RS2 


this work 


0.028 


0.672 


1.329 


0.345 


-0.015 






Ref. Jl7| 


0.030(3) 


0.678(10) 


1.330(17) 


0.349(5) 


-0.034(30) 


XY 


RSI 


this wojjk 


0.029 


0.667 


1.318 


0.343 


-0.001 






Ref. Jiq 


0.034(3) 


0.669(1) 


1.316(1) 


0.346(1) 


-0.007(6) 


Heisenberg 


RSI 


this work 


0.028 


0.697 


1.379 


0.369 


-0.092 






Ref. 


0.034(3) 


0.705(1) 


1.387(1) 


0.364(1) 


-0.115(9) 
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